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|"t ! ' We study the Higgs sector in the Next-to-Minimal Supersymmetric Standard 

Model with and without explicit CP violation, focusing on the case of weak scale 

£NJ \ expectation value of the singlet field. We scan a wide range of the parameter space 

to find out allowed regions by requiring that the electroweak vacuum be the global 
minimum of the effective potential and that the neutral Higgs bosons with moderate 
gauge coupling be heavier than the lower bound on the Higgs boson in the standard 
| model. Among the allowed parameters, some sets admit the situation in which the 

■ light Higgs bosons couple with the Z boson too weak to be found in present collider 
. experiments. For such parameter sets, we find an upper bound on the charged Higgs 
\ mass which is reachable in LHC. 

^ : 

^ 1 Introduction 

The search for the Higgs boson is one of the most important issue of high-energy 
particle physics because the Higgs boson is the only unobserved particle in the 
Minimal Standard Model (MSM). The results of the LEP 2 experiment place the 
limit on the MSM-Higgs mass: m h > 114.4GeV at 95% CL [1,2]. Although there 

■ are some theoretical restrictions to the Higggs mass, it cannot be predicted in the 
MSM framework, since the Higgs self-coupling is a free parameter. Supersymmetric 
extensions of the MSM, which are motivated to solve the hierarchy problem, limit the 
range of the Higgs mass because of the self-coupling given by the gauge couplings. 
Among such extensions, the Minimal Supersymmetric Standard Model (MSSM) 
has been well studied and is known to give the upper mass bound of the lightest 
Higgs boson, m^ < mz at the tree level. This bound seems somewhat severe but 
it is modified by radiative corrections, shifting it to about mh < 135GeV at the 
two- loop level, mainly by the loops of the top quark and squark [3]. 

The MSSM contains a /^-parameter in the superpotential. It enters the Higgs 
potential with the soft scalar masses to determine the vacuum expectation value 
(VEV) of the Higgs fields. Then \x must take a value of weak scale, which is much 
smaller than the GUT scale or Planck scale. There is no a priori reason for fi to 
have such a small value. One solution of this so-called /U-problem is to substitute 
a VEV of an extra gauge-singlet field for the parameter fi. The NMSSM is among 
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such models which have a gauge-singlet Higgs superfield N. The superpotential of 
the model contains 

W = -\NH d H u - ^kN 3 , (1) 

in addition to the MSSM terms with \i = [4]. We adopt the ^-symmetric version 
of the superpotential so that it does not contain any dimensional coupling. The // 
parameter of the MSSM is generated as \i = A (N) . The NMSSM behaves like the 
MSSM in the limit of (N) 3> v with \{N) and k(N) fixed, for which the singlet 
decouples. Here v = \J {H^) 2 + {H u ) 2 is the VEV of the Higgs doublets. Since 
the superpotential has no dimensional parameters, the VEVs of the Higgs fields 
are determined by the mass parameters in the soft-supersymmetry-breaking terms 
together with the couplings. If all the mass parameters in the Higgs potential are of 
weak scale and all the couplings have moderate values about 0.1 — 1, it is natural for 
the singlet to acquire a VEV of the same order as v. Then we expect new features 
which are absent from the MSSM, and our main concern is in the Higgs sector of 
this case. 

The NMSSM contains three CP-even neutral Higgs bosons (Si, S2, S3), two CP- 
odd bosons (Ai,A%) and a pair of charged bosons in the CP-conserving case. 
The spectrum of the CP-conserving model has been studied in [5-7]. In contrast 
to the MSSM, the three CP-even scalars can mix up to form mass eigenstates of 
small mass with very small gauge coupling [7]. Such a light Higgs boson cannot be 
produced at lepton colliders so that it is not excluded even if its mass is smaller 
than 114GeV. We shall refer to this situation as light Higgs scenario. As we see 
below, such a light Higgs situation is realized for weak scale (N) and small k. A 
similar situation has been observed in the MSSM, when a large mixing among CP 
eigenstates is caused by the CP-violation in the squark sector, which is characterized 
by the imaginary part of the product of \i and the A-term [8]. Then the gauge 
coupling of the lightest scalar is so small that it can escape the lower bound on 
the Higgs mass [9]. While the same mixing is also expected in the NMSSM, it 
contains another source of CP violation in the tree-level Higgs sector. Such a CP 
violation has been studied in several cases where it is caused spontaneously [10] and 
explicitly [11] in some special situations. We give the one-loop formulation with 
all possible CP phases including squark sector in a manner independent of phase 
convention. In this formulation, one can easily arrange the phases in such a way that 
the phase relevant to the neutron Electric Dipole Moment (nEDM) is suppressed 
while retaining those which affect the mixing of the Higgs bosons. We investigated 
the mass spectrum and couplings in the presence of such a CP phase. 

Another aspect of the NMSSM Higgs sector is that the Higgs potential contains 
cubic terms including the singlet field. Although these terms must be constrained 
not to generate the global minimum of the effective potential different from the 
electroweak vacuum, they are expected to make stronger the first-order phase tran- 
sition at high temperature. In this sense, this model is more suited for electroweak 
baryogenesis than the MSSM, which requires a light stop with mass less than the 
top quark mass for the strongly first-order phase transition [12]. We pointed out 
that the CP violation in the squark sector weakens the electroweak phase transition 
(EWPT) caused by a light stop in the MSSM [13]. We expect that, in contrast to 
the MSSM, the CP violation in the Higgs sector of the NMSSM will not weakens 
the EWPT, while supplying sufficient CP violation to generate chiral charge flux, 
which is the source of baryon asymmetry. The effects of this CP violation on the 
phase transition will be discussed in our forth-coming paper. 
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This paper is organized as follows. In Section 2, we analyze the NMSSM Higgs 
sector at the tree level and explain how to obtain restrictions on the parameters in 
the model. There we derive the upper and lower bounds on the mass of the charged 
Higgs, which are trivial in the MSSM limit but important in the case of weak scale 
(N). In Section 3, we show the one-loop formulas for the mass-squared matrix of the 
neutral Higgs bosons and a mass of the charged Higgs boson. Section 4 is devoted 
to the numerical results for the parameter search. The spectrum condition divides 
the allowed parameter sets into two classes. The first one is the MSSM-like allowed 
parameter sets where all the Higgs bosons are heavier than the 114 GeV, and the 
second one corresponds to the light Higgs scenario. The study of the CP violation 
is described in Section 5. The formulas used to define the effective potential and to 
calculate the mass matrix are summarized in Appendices. 

2 Tree-level Higgs sector 
2.1 Higgs potential 

In this Section we analyze the tree-level Higgs sector. The NMSSM has the super- 
potential with the singlet superfield iV [7], 

W = f d H d QD c - f u H u QU c - XNH d H u - ^N 3 , (2) 

where Q, D c and U c denote chiral superfields containing quarks, and H d and H u 
contain the Higgs doublets required in the MSSM. The couplings A and k for the 
singlet are in general complex numbers. We consider ^-symmetric version of the 
superpotential, so it does not contain any dimensional coupling. 

In addition to the supersymmetric Lagrangian, the low-energy NMSSM contains 
the soft-SUSY-breaking terms, 

Asoft = - m\§\<S> d - ml<S>i$ u - m 2 n n*n - m\q\q L - rn^d ] B d R - m\u R u R 

- \[f d A d )$ d q L d* R - UuA u )$> u q L u R - \A x n$ d <f> u + h.c.} 

- {m' n n 2 + ^kA k u 3 + h.c), (3) 

where q~L, d R and u R are the squark fields and <fr d , an d n are the Higgs fields. 
Although the n 2 -term breaks the global ^-symmetry, the n 2 -term does not exist 
in the simple supergravity model. Hence we shall not include the n 2 -term in the 
following. We assume that all the dimensional parameters in £ so ft have their values 
near weak scale. 

In the limit of k = 0, the global ^-symmetry is elevated to U(l) PQ symmetry 
(n — > e ia n and $ u — > e~ ia & u ). Then the pseudoscalar component of the singlet 
becomes a Nambu-Goldstone (NG) mode when the singlet acquires a VEV. At 
small k, the PQ symmetry is slightly broken and relatively light axion is expected 
[14]. The spontaneous breakdown of the global ^-symmetry causes the domain 
wall problem [15]. If we introduced ^-breaking linear or bilinear terms into the 
superpotential, this problem could be solved [16]. It inevitably spoils the nature 
of the superpotential with no dimensional parameter as the countervalue. So we 
assume that this symmetry is broken by some higher dimensional operator which 
becomes irrelevant in low energy physics, at the early stage of the universe far before 
the EWPT. 
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The tree- level Higgs potential is composed of three parts, V = Vp + Vd + V^ D f t ; 
V F = \Xn\ 2 (^d + + \eijXmi + ku 2 \ 2 , (4) 

Vd = ^±^(*J*d - *I*«) 2 + f (*S*«)(*l,*d), (5) 
y soft = m?$t$ d + ?n^t $ u + m 2 N \n\ 2 - {e^XA^^l + ^kA k u 3 + h.c.). (6) 

Here we expand this potential around the VEVs which are represented by v d , v u , 
v n and the phases and (p. The parametrization of the scalar fields are as follows, 



(vd + h d + ia d ) 



n = -^=e lv> (v n + h n + ia n ). (8) 

The condition for the scalar potential to have an extremum at the vacuum is 
that the first derivatives with respect to the Higgs fields evaluated at the vacuum 
vanish: 

1 /dV \ „ 2 V u V n gl+g(, 2 2n . I A I% 2 . 2n . ^ v u^l fa . 



1 / dV \ „ 2 V d V n g 2 +g 2 f 2 2n . 2 , 2\ ■ nv dvl 

- m 2 - R x — (v d - vj + —(v d + vj + 



v u \oh u v u 8 2 2 w u 



(10) 



i / \ _ _ 2 _ _ R . w _ p ^ |A| 
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=— ( air ) = ™N ~ R\— ~ RnVn + hrivj + <) + |«|X + ^«^«, (U) 

v n \oh N / v n 2 

i /dv \ i /avb\ r i 2 ^ 

=— ( fl— > = — ( fl— > = hv n - 7> n> ( 12 ) 

u« \da d / v d \da u / 2 

0= — — )=- f A \- I K v n +lv d v u , (13) 

u n \da N / f n 

where 

U = Re[\K*e^ e - 2 ^}, 1 = Im[Arc*e^- 2 ^], (14) 

R x = -^=Re[XA x e^}, I x = -±=Im[\A x e^}, (15) 

R K = -^ReM K e <3 % /« = -^Im[d K e ,3 % (16) 

and (• • • ) denotes the value evaluated at the vacuum. These conditions are called 
tadpole conditions in the sense that the conditions make the tadpole diagrams van- 
ish, if we set the Higgs fields to their VEVs. Note that 1Z and I are dimensionless 
parameters and R x , R K , I x and I K have dimension one. The phases only appear 
through the three combinations (14)-(16) given above. Hence our formulation so 
far does not depend on the convention of the phases. From (12) and (13), we obtain 
two conditions 

= ~ Iv n, Ik = • (17) 

z z v n 

Because of the two tadpole conditions on the three CP violating parameters X, 
I x and I K , only one of them is the physical one. When we introduce complex 
parameters, we have to manage them to satisfy tadpole conditions (17). 
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2.2 The mass and couplings of the Higgs scalars 

The mass matrix of the neutral Higgs scalars, which is defined by the second deriva- 
tive of the Higgs potential evaluated at the vacuum, has the following structure, 



M' 



M 
(M 



2 

S 

2 \T 
SP) 



Ml 



(18) 



where the basis is (h T a T ) 
M 2 are given by 



(hd h u h n ad a u a n ). Here the block components of 



/ R\v n tan + m 2 z cos 2 
— TjTZv 2 , tan 



Mi 



-R\v n — m 2 z sin cos 
+\Hvl + \\\ 2 v d v u 



-R\v n — m 2 z sin cos R\v n cot + m 2 z sin 2 



M 2 P 



+\Kv 2 n + \\\ 2 v d v u 

-R\v u + \ \\ 2 v n v d 
+Kv n v u 

(Rx - Uv n /2)v n tanf3 
(Rx - Kv n /2)v n 
(Rx + Tlv 



-7}Kv 2 cot/? 

-R\Vd + \A 2v nV u 

+TZv n v d 



-R x v u + \X\ 2 v n v d \ 
+lZv n v u 

-Rxv d + \X\ 2 v n v u 
+Tlv n Vd 



(Rx - Tlv n /2)v n 
(Rx - Kv n /2)v n cot 
(Rx +TZv n )v d 



M 2 



SP 







\lv n v d 



RkVu 

+2\k\ 2 v 2 h 

(19) 

(Rx + Tlv 

n )Vu 

(R X + TLVn)Vd 

(20) 

%Tv n v u 

(21) 



2lv n Vd 

-2Iv d v u 



where we use the tadpole conditions (9)-(ll) and (17) to express the scalar soft 
masses and Ix, I K i n terms of the other parameters. In the following, we also adopt 
the usual conventions tan/3 = v u /vd and v 2 = v^ + v 2 . It is worth emphasizing that, 
if CP violation at the tree level is turned off (i.e. 1 = 0), the scalar mass matrix 
becomes block diagonal. 

The three pseudoscalars contain one NG mode, which is isolated by the (5 rota- 
tion 

cos (5 sin P 
- sin cos 
1 



a 



U(0) 



G 
a' 









(22) 







where G is the NG mode which would be eaten up by the gauge bosons. After 
isolating the NG mode, the mass term of the neutral Higgs bosons is given by 



1 



( fc T 



a 



iT 



M 2 S 



(M 



SP 



My 

M 2 p 



h 

a' 



(23) 



where 



Ml' 



M 2 SP ' 



(2R X - TZv r , 



Un > sin 2/3 

(R x + Kv n )v 



| sin ,5 
| cos 
-sin 20 



2v n v, 



(Rx+7lv n )v 
% f- sin 20 + 3R K v n - Kv 2 sin 20 



(24) 
(25) 



From the equations (24) and (25), we are convinced that one of the pseudoscalar 
becomes massless in the limit k = (X = 1Z = R K = 0). This holds even when 
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we include radiative corrections. We define Ai' 2 as the mass matrix of the neutral 
Higgs bosons after extracting the NG mode. Then the mass of the neutral Higgs 
bosons are obtained by diagonalizing the mass matrix Ai' 2 by orthogonal rotation 
T Ai' 2 = diag(m| i m 2 ^ m 2 ^ m 2 ^ m^), where we define the matrix O in such a 
way that the eigenvalues satisfy 

™hi < m l 2 < m l 3 < m li < m t 5 - ( 26 ) 

Without CP violation (X = 0), the mass eigenstates are also CP eigenstates Si and 
Aj, where Ai has vanishing gvvh coupling (V represents the W and Z bosons). 

Similarly, the charged Higgs mass m H ± is obtained by the (3 rotation of the 
charged Higgs mass matrix, 

2 1 / d 2 V 



sin /3 cos /3 \ d(f)t dcj) u 



m 2 W ~ 7^1 V + (2i?A - Kvn)-^- (27) 



where 



We use this equation in order to substitute mjj± for R\ by 

„ 1 . 9 sin 2/3 1 . . 

R\ = ~™ + ~Kv n , 28 

2 v n 2 

m 2 = m 2 H± -m 2 w + ^\X\ 2 v 2 . (29) 
For example, the CP-odd components of the mass matrix is written as 

V \rh 2 ^s\n2(3 + \llv n v \m 2 {^ sin2/3) 2 - \Kv 2 sin2/3 + 3K K v n 

(30) 

Now we have seven mass eigenstates, but we have no obvious bounds on the 
eigenvalues, like the upper (lower) bound on the lightest (heaviest) Higgs scalar 
in the MSSM at the tree level. Instead, without CP violation, the inequality 
det (to 2 — .Mp) < implies that m\ < rh 2 < vnr A . It is difficult to derive such in- 
equality for the CP-even scalars, but in the limit of m 2 3> Vq, f 2 , we have an approx- 
imate relation det (m 2 — -M|) < 0, which implies that < m| 2 < rh 2 < m 2 S3 . 
These relations explain the pattern of the mass eigenvalues in the case of heavy 
charged Higgs boson, with the help of Tr.A/f | and Tr.Mp, which constrain the sum 
of the masses. 

Although the singlet Higgs fields h n and a n do not couple to the gauge bosons, 
all the mass eigenstates of the neutral Higgs boson can interact with the W, Z 
bosons and fermions, because the singlet fields are mixed with the doublet fields. 
The couplings of the charged Higgs boson with the gauge bosons and quarks are 
identical to those in the MSSM. In particular, the VVh-, Zhh- and 66/i-vertices 
are important for the study of Higgs production and decay events in colliders [2]. 
At the LEP-type e + e~ collider, the dominant production processes of the neutral 
Higgs bosons are the Higgs strahlung process associating with Z boson, the W- 
fusion production, and the pair production processes if +mh 2 < y/s- As for the 
decay of a light Higgs boson, h — > bb is the main mode, with subleading h — > t + t~ 
mode. In both cases, the correction factors to the relevant Yukawa couplings are 
the same with each other which characterize deviation from the MSM. In addition 
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to the processes, the gluon fusion and the Yukawa processes become important at 
high-energy hadron colliders [17]. 

It is straightforward to read off those vertices from the kinetic terms of the Higgs 
bosons, and from the Yukawa coupling terms: 

Cvvh = 92mw9vvh i (w+W-» + 2c ^2q w Z » Z ^) ( 31 ) 

£zhh = 92 a (hdd/iad - h u d^a u )Z^ (32) 
2 cos t>w 

gzhihiZ^ihid^hj), 



2 cos 9w 

tbbh = -^^^dbbhi + rfdwhjbhi, (33) 
where the correction factors to the couplings are 

QVVhi = On cos + 2i sin /3, (34) 
Qzhihi = o {{0 4i 2 j - Oij0 2 i) cos/3 - (OuOij - OijOu) sin/3} , (35) 

9bbh t = On— — r, g^ hh . = -O u tan p. (36) 

Then the mass eigenstates and the gauge eigenstates are related as = Ouhi, h u = 
2 ihi and a = Ouhi, where i is summed from 1 to 5. The Zh{hj coupling vanishes 
for i = j, because of the antisymmetric derivative nature of its definition. This 
coupling also vanishes when both the Higgs bosons are either scalars = 1,2,3) 
or pseudoscalars = 4,5). However, all the gzhih 3 are expected to have nonzero 
values in the CP-violating case, because of the mixing of these CP eigenstates. The 
equations (34)-(36) have the same structure as those in the MSSM [8], except for 
the mixing including the singlet. 



2.3 Constraints on the parameters 

The NMSSM has more parameters than the MSSM. Our main concern is to search 
for allowed parameters in the case of weak scale v n , for which we expect new features 
in the spectrum and coupling of the Higgs bosons, as well as the phase transitions 
at finite temperature. In order to select the allowed parameters, we impose the 
following two conditions on the model: 

(1) the vacuum condition, which requires that prescribed vacuum is the global 
minimum of the effective potential. This also requires that all the masses- 
squared of the scalars including the sfermions be positive. 

(2) the spectrum condition, which requires that the mass of the Higgs boson with 
its couplings to the vector boson |<7yv/i| larger than 0.1 be heavier than the 
bound 114GeV. 

The constraint on the gauge coupling is the most stringent so that we examine the 
other couplings for the allowed parameters later. Since the mass matrix of the Higgs 
bosons receives large radiative corrections, we need some numerical studies to figure 
out the results of the spectrum condition, which will be presented in Section 4. Here 
we attempt to find analytic form of constraints obtained from the vacuum condition 
at the tree level. 
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In the MSSM, the global minimum of the tree-level potential is always located 
at the vacuum, as long as the D-flat direction is raised by the soft terms and 
the tadpole conditions are satisfied. Although the Higgs potential in the NMSSM 
is bounded from below by the .F-terms, the prescribed vacuum is not always the 
global minimum of the potential, even when the tadpole conditions are satisfied. 
This is because the trilinear terms in £ so ft, which give negative contributions to the 
potential, make some point different from the vacuum to be the global minimum. 
We must exclude such a parameter set which yields an unwanted global minimum. 
A necessary condition for the correct vacuum is that the mass-squared of all the 
scalars be positive. In the CP-conserving case, this applied to the pseudoscalars 
implies that det M.' p 2 > 0, hence, 

rh 2 (--Kv 2 sin 2/3 + R R v n ) > -Tl 2 v 2 v 2 n . (37) 

This requires that each factor in the left-hand side has the same sign, and gives the 
lower bound on the charged scalar mass, which is large enough for m 2 > 0. 

Another necessary condition is that the value of the scalar potential at prescribed 
vacuum be lower than that at the origin. Now at the tree level, V(0) is zero and 
that at the vacuum is 

t r I 1|\|2 2 2 1| |2 4 1 2 2 2 r,/3 ^ 2 2 • 2 tf 

V \ vacuum = ~ 4 1*1 V n V ~ 4M V n ~ ^ Z V COS 2(3 - -m w V Sill 2(3 

+ \m 2 H± v 2 sin 2 2/3 - \lZv 2 n v 2 sin 2/3 + -R K v z n , (38) 
8 8 o 

where we use the tadpole conditions (9)-(ll) to eliminate the soft masses of the 
Higgs fields and use m H ± instead of R\. Then requiring that V| vacuum < gives 
m H ± the upper bound 



m 2 H± < 2\X\ 2 v 2 -^— + 2\k\ 2 ^-^— + m| cot 2 2/3 + m 2 w 



A v 3 1 
—R — 

1 sin 2/3 3 K v 2 sin 2 2/3 



+ n <— "^5 ~ oRk^-tt^Z- ( 39 ) 



This bound is available regardless of whether there is CP violation or not. The 
charged Higgs mass mjj± is not constrained in the MSSM-limit where \v n and Kv n 
fixed for v n — ► oo [6], because of the infinitely large term. So this condition 
becomes important in the pure-NMSSM parameter set, i.e., when v n is not so large 
and A and k are not so small. Fig. 1 show the tree-level charged Higgs mass bounds 
for an example of pure-NMSSM parameter set. The solid line which shows the 
upper bound suggests that the charged Higgs boson must be lighter than 400 GeV. 
Hence the pure-NMSSM parameter set predicts the charged Higgs boson with the 
mass accessible to the LHC [18]. The dashed line indicates the lower bound on the 
charged Higgs boson. It ends near A K = because the left hand side of (37) becomes 
negative. The consistency of the model sets both the upper and lower bounds on 
A K , which are of weak scale for a weak scale v n . In particular, the same sign of k 
as A K is favored. 

Once we fix m#± and A K , the condition (39) excludes the elliptic region in the 
(A, «)-plane, whose area vanishes in the MSSM-limit. As we see below, this excludes 
a large portion of the parameter space for the case of weak scale v n . 
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m H L [GeV] ""si [GcV] 
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A K [GeV] 4. [GeV] 

Figure 1: The mass bounds on the tree-level charged Higgs boson as a function of A K for 
tan/? = 5, v n = 300GeV, A = 0.1 and k = —0.3 (left-hand plot) and k = 0.3 (right-hand 
plot), respectively. The solid line shows the upper mass bound of charged Higgs boson 
and dashed one shows a lower bound. 



3 One-loop Effective Potential 

We now evaluate the one-loop contributions to the tadpole conditions and the spec- 
trum of the Higgs bosons. As we mentioned in the previous Section, the basic idea 
is almost the same as in the tree-level analysis, so we show the results briefly We 
analyze the effective potential of the Higgs fields by taking into account the one- 
loop contributions from the gauge bosons and the third generation of the quarks and 
squarks. The corrections from the leptons, the other quarks and squarks can be ig- 
nored, because of their small Yukawa couplings to the Higgs fields. The corrections 
to the Higgs potential are given by 



= A q V + AsV + A„V, 



(40) 



where 



A q V 



A g V 



N c 
16vr 2 



q=t,b 



m 



2^2/ 1 m q 

V og w 



N c 

327T 2 

3 

64vr 2 



E E (™ 



q=t,bj=l,2 
(ml? 



1} 



rrir 



log 



M 2 



l0g JP~2 



2 (fh 



w 



log 



m 



w 



M 2 



(41) 
(42) 
(43) 



The field-dependent masses m\ are listed in Appendix A. The tadpole conditions 
and the mass matrix are obtained by calculating the derivatives of the effective 
potential at the vacuum with these corrections. In particular, the imaginary parts 
of the tadpole conditions are 



h + A A 



1 



-lv n , I K 



2 V n 



(44) 



9 



where I K does not receive the loop correction at this level. The correction AI\ 
comes from the squark loops 

N c 



Ah 



where l„ is defined by 



16vr s 



Yl \y<l\ 2l Q [/( m L m |) - !] "n, 



q=t,b 



and 



(45) 



f(ml,mf) 



Am 2 



2/1 m i 
m i ( lo S 



m ( lo 



777-2 

F 



(46) 



with Am 2 = m\ — 

Here we show the corrections to the mass matrix in detail. The couplings to 
the gauge bosons and the quarks receive the loop effects, through the orthogonal 
matrix O which is determined from the corrections to the neutral mass matrix M. 
The scalar part M 2 S of the neutral mass matrix receives the loop contributions, 

AM% = A q M 2 s + A g Ml + A q -M% (47) 

from the quarks, gauge bosons and squarks. The quark loops and gauge loops 
contribute to the upper-left 2x2 elements only, 



m b 
IP 



1 1 9 9 1 
\yb\ ml log- 






\yt\ 2 m 2 log j 




(48) 



A 9 M 2 s 



32vr 2 



rrir 



m z log ^ + 2m w log 




cos 2 (5 cos j3 sin f3 
cos /? sin (3 sin 2 /3 




(49) 

because they do not couple to the singlet in the tree-level potential. The squark-loop 
contributions are 



N c 
16vr 2 



T S f(ml,ml)+Y: 



3=1,2 



i=i,2 



'9m 2 
9/i 



^ 

5/i 



' 9m 2 



' dm- 
h ± 

dh 



l0g # 



T 



m~ 



(50) 



(51) 



where the matrices T, B and the list of derivatives of the field dependent masses, 
are given in Appendix B. 

M 2 P and M 2 SP are not affected by quark and gauge loops but receive the con- 
tributions from squark loops 



AMI = &<}M 2 P , AM 2 SP = AqM 2 SP , 



where squark-loop contributions are 

N c 



A-M 2 P 



A~ h M 2 P 



16vr 2 

N c 
16vr 2 



3=1,2 



W ' \ Q a 

3=1,2 




(52) 

(53) 
(54) 
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and 



2 

SP 



A~ b M 2 SP 



N c 
16vr 2 

N c 
16vr 2 



-r 



SP 



771 



77? 



6 2 



+ E 

.7=1,2 



+ E 

.7=1,2 



' dm- 
dh 



'dm 2 
b _i 

dh 



' dm 2 
da 

I dm- 
h 

da 



log 



2 

m~ 
l j_ 

M 2 



m~ 



log 



M 2 



(55) 



(56) 



The NG mode can be extracted from M/p and M SP by the same orthogonal trans- 
formation as at the tree level. The loop contributions to the off-diagonal matrix 
M%p are proportional to the CP-violating parameters in the squark sector, I q which 
is defined in (68). 

The charged Higgs mass is not affected by the singlet field but receives the 
contributions from the gauge, quark and squark loops, 



Am 2 H± 



A g m 2 H ± + A q m 2 H± + A q m 2 H ±, 



(57) 



where the detailed form of the each terms are given in Appendix C of our previous 
paper [13], except that /j in the MSSM is replaced by Av n e lip /v2. R\, which is 
determined from equation (28), receives loop corrections through the corrections to 
the charged Higgs mass. 



4 Parameter search 

In this Section, we search for the allowed parameter region numerically with the 
one-loop corrections in the CP-conserving case. The allowed parameter sets are 
selected by requiring the two conditions discussed in Section 2.3. For this purpose, 
we scanned the parameters in the Higgs sector within the following region: tan (3 = 
3 - 20, v n = 100 - lOOOGeV, 100 < m H ± < 5000 GeV and -1000 < A K < 
GeV. Since we can always make A positive without loss of generality, (A, K)-plane 
is scanned for < A < 1 and — 1 < k < 1. As for the squark sector, we adopt a 
small value of the ^4-term, At = Ab = 20GeV, for the squark fields not to acquire 
nonzero VEV, and three patterns of the soft masses: the heavy-squark scenario 
with (rriq, mi R ) = (lOOOGeV, 800GeV), the light-squark-1 with (m,q, m^ R ) = 
(lOOOGeV, lOGeV) and the light-squark-2 with (m q , mQ = (500GeV, lOGeV), 
where m q (mf; R ) denotes the doublet (singlet) soft mass. For simplicity, we put 
m b R = m t R - First of all, we pick up a set of parameters except for A and k and 
exclude the regions in (A, «;)-plane where the effective potential at the origin is 
smaller than that at the vacuum and the spectrum condition is not satisfied. Within 
the remaining region, we perform the numerical search for the global minimum of 
the effective potential and exclude the points for which the minimum is different 
from the vacuum. 

In Fig. 2, we present a typical example of the allowed parameter region for 
tan (3 = 3, v n = 200 GeV, m H ± = 400 GeV and A K = -200 GeV for heavy-squark 
scenario. The white region indicates the allowed parameter region in (A, K)-plane. 
Although negative k is favored for negative A K as explained in Section 2.3, we 
study the range of k as —1 to 1 since it is not excluded completely. Within the 
dark gray elliptic region, the effective potential at the origin is smaller than that at 
the prescribed vacuum. This is predicted by (39) which was derived from the tree- 
level potential. The broad light gray region is excluded by the spectrum condition. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



A 

Figure 2: Allowed parameter region as a function of A and k for tan j3 — 3, v n = 200 GeV, 
rriH± = 400 GeV and A K = —200 GeV for heavy-squark scenario. Allowed parameter 
region are shown in white. 



Within the upper allowed region near k = 0, the light Higgs scenario is realized, 
while the lower allowed region yields the lightest Higgs boson heavier than 114 
GeV. The narrow black region right to the elliptic one is excluded, since the global 
minimum is located at a point different from the prescribed vacuum. For these 
excluded parameters, the global minimum is located at v = and a large v n , which 

9 i 1 2 

depends on m N , R K and |«| . 

We explain the dependence of the allowed regions on the parameters. In Fig. 3, 
we show the same as Fig. 2, but with (a) v n = 500 GeV, (b) A K = 500 GeV, (c) 
light-squark-2 scenario and (d) tan (3 = 5 and m jj± = 600 GeV. As shown in Fig. 3 
(a), for a larger v n , the allowed region of the light Higgs shrinks and the allowed 
region of the heavy Higgs spreads out to small A value. The allowed parameter set 
with small A corresponds to the MSSM limit. Although not shown in the graph, the 
similar behavior is observed, when the charged Higgs boson becomes heavier. For 
larger A K , the allowed region becomes smaller. It is shown in Fig. 3 (b) and was also 
expected from Fig. 1. If we choose small m« and rriq R , Fig. 3 (c) and the same plot 
for the light-squark-2 scenario show weak dependence on the squark soft masses. 
Fig. 3 (d) indicates that for larger tan/3, the allowed region with light Higgs gets 
narrower. At tan (3 = 20, the allowed region becomes point like at k = 0, where one 
of the pseudoscalar is always massless. 

Now, we consider the detail of the Higgs mass spectrum and couplings. In Fig. 4 
and 5, we show the behaviors of the Higgs masses (Fig. 4 left) and the couplings 
of the three lightest Higgs bosons to the massive gauge bosons (Fig. 4 right), to 
the bottom quarks (Fig. 5 left), and the Z/i/i-couplings (Fig. 5 right), as a function 
of k for the same parameters as Fig. 2 with A = 0.9. These correction factors to 
the coupling constants are defined in (34)-(36). As seen from Fig. 2, the allowed 
region along A = 0.9 is divided into two parts, one of which corresponds to the light 
Higgs scenario with small \n\. The range of k satisfying the spectrum condition 
can be read from Fig. 4 to be —0.33 < k < —0.05, some part of which is excluded 
by the vacuum condition. Within the other allowed region for —1.0 < k < —0.8, 
all the Higgs bosons are heavier than 114 GeV. We refer to this parameter set as 
the heavy Higgs scenario. In the scenario, the relatively heavy Higgs boson /13 is 
almost decouple from the theory because its couplings gvvh 3 and gbbh 3 are small. 
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Figure 3: The same as Fig. 2, but with (a) v n = 500 GeV, (b) A K =500 GeV, (c) light- 



squark-2 scenario (i.e. rriq = 500 GeV, mi 
and m H ± = 600 GeV. 

[GeV] 
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K 



10 and A t = 20 GeV) and (d) tan/3 = 5 



9vvh 3 
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Figure 4: The Higgs masses (left) and the couplings (right) to the V-bosons as a function 
of k for the same parameters as Fig. 2 with A = 0.9. In the left plot, solid line, broken 
line and dotted line represent m^, m^ and m ha , respectively. In the right plot, solid line, 
dashed line and dotted line represent gvvhv 9vvh 2 an< ^ dvvh^i respectively. 
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Figure 5: The couplings of the three lightest Higgs bosons to the bottom quarks and the 
Z boson as functions of k for the same parameters as Fig. 4. In the left plot, solid line, 
broken line and dotted line represent g 2 hhhx = (g§ bhl ) 2 : + {g bhhl f ', g bbh2 , gl bh3 , respectively. In 
the right plot, solid line, broken line and dotted line represent g 2 Zhlh2 , 9zhih 3 an d 9zh 2 h 3 i 
respectively. 



The lighter two Higgs bosons are both CP-even scalars, hence in the scenario, the 
NMSSM behaves like the MSSM. Particularly the lightest Higgs boson is heavier 
than 120GeV and its coupling to the bottom quark is not so large, g bbh < 2.2, which 
escapes the experimental bounds as the MSSM does. It is difficult to distinguish 
the NMSSM from the MSSM in the heavy Higgs scenario. 

In the light Higgs scenario, the light Higgs bosons cannot be observed in the 
collider experiments. The main processes for producing the Higgs bosons at the 
LEP 2 are the Higgs strahlung process and the VF-fusion process, but both processes 
include small FF/i-couplings (shown in the right plot of Fig. 4) and hardly produce 
the light Higgs bosons. Though the pair production becomes important when the 
total mass of the lightest scalar and the lightest pseudoscalar bosons is under the 
LEP 2 threshold, they cannot be created in pair [19] because of so small coupling 
9 Zh\hz (9zh 1 h 2 ^ shown in the right plot of Fig. 5. The Yukawa processes can be 
dominant where these processes are suppressed but the Yukawa processes cannot 
be observed unless the bbh-coup lings are fairly large [20]. In the left plot of Fig. 5, 
the Yukawa couplings g bb h are almost unity in the light Higgs region. With these 
couplings the light Higgs bosons h\ and hi produced in the Yukawa processes are 
hardly observed. These arguments are essentially the same even if we consider the 
hadron colliders. The SM-Higgs search experiments at the Tevatron mainly consider 
the strahlung process with W(Z) in the low mass range 110 < nth < 140 GeV, and 
the results show that the signal efficiency is too small to find out the Higgs boson 
[21]. At the LHC, the gluon fusion and the IF-fusion processes become important. 
Though the h\ and hi production in IF-fusion process is suppressed by their small 
gauge couplings, their production in gluon fusion process is not suppressed. However 
it is in general difficult to select the gluon fusion events from the backgrounds, 
because the process produces the Higgs boson only. Therefore /13 is expected to be 
first detected in the collider experiment. Then the model is very similar to the SM 
with the Higgs mass which can be as heavy as 150GeV. 
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5 Effects of the CP violation 



Now we turn on CP violation and study its effects on the spectrum and gauge 
couplings of the Higgs bosons. Such an analysis in the MSSM by considering the 
CP violation in the squark sector was done in [8]. Here we focus on the tree- 
level CP violation which does not exist in the MSSM and is characterized by the 
parameter X defined in (14). In addition to Z, the complex parameters in the Higgs 
sector entering the mass-squared matrix are summarized in 1Z, R\ and R K . When 
all the parameters are real as in the previous Section, one can freely assign their 
values. However, if some of them are complex, we must arrange the parameters and 
the phases and ip to satisfy the tadpole conditions (44). Before presenting the 
numerical results, we explain how we parameterize the CP violation. 

First of all, R\ is fixed by the charged Higgs mass m H ± from the equations (28) 
and (57). The remainders, TZ, R K and I, are determined by A, k, A k , 8 and ip. 
Here A, n and A K are complex numbers and some of the phases are redundant. We 
denote the phase of A, k and A K as 5\, 5 K and 8a k , respectively. Among those, 
independent phases are corrected into 8' K = 5 K + 3tp, 5a k and <5edm = + + ^P, 
which is effective to nEDM if the gaugino masses and A q are real. The counterpart 
of the <5edm in the MSSM is the phase of the /x-term plus 6. Suppose that we first 
give |A|, \k\, <5edm and 5' K , from which 1Z and I are determined. Since I K is fixed by 
(44), one can take any absolute value of A K but not its phase. In particular, R K is 
given, without specifying 5a k , by the relation 

where we use the second equation of (44). Not that, the \nA K \ must be larger than 
3IvdV u /\f2v n in order for R K to have a value. In the following, we adopt the phase 
convention so that the relevant phase to the nEDM vanishes. Therefore, 5' K is the 
only phase that we can freely assign its value. 

Finally we show the effects of S' K on the masses and couplings of the Higgs 
mass eigenstates with <5edm = 0. For illustration, we take two parameter sets with 
k = —0.2 and k = —0.9, while the other parameters are the same as in Fig. 4, and 
plot ^-dependences of the masses and couplings in Fig. 6 and Fig. 7, respectively. 

In order to compare with the former result, positive R K in the equation (58) was 
chosen. As shown in Fig. 6 and Fig. 7, the phase dependences in the light Higgs 
scenario (k = —0.2) are smaller than in the heavy Higgs scenario (k = —0.9), since 
\k\ is smaller for the light Higgs scenario. In Fig. 6, the next-to-lightest Higgs 
boson /12 becomes lighter than 114 GeV for S' K /ir < 0.31, hence the large 5' k is not 
allowed by the spectrum condition. In Fig. 7, the model is excluded for 5' K /ir > 0.24, 
because of the small Higgs mass with moderate gauge coupling and also excluded 
for S' k /tt > 0.39 where the prescribed vacuum becomes unstable. We could not find 
an example in which an excluded parameter set is converted to be allowed as a 
light Higgs by introducing the CP violation in the tree-level Higgs sector. However 
the light Higgs boson is realized even for the heavy Higgs parameter sets when we 
introduce CP violation in the squark sector, in the same way as in the MSSM. 

6 Summary 

We investigated the spectrum and coupling constants of the Higgs bosons in the 
NMSSM, when the vacuum expectation value of the singlet is of the weak scale for 
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igure 6: The same as Fig. 4 but with k = —0.2 (light Higgs scenario) and Sedm = 
function of 5' . 
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Figure 7: The same as Fig. 6 but with k = —0.9 (heavy Higgs scenario). 



a wide range of all the parameters in the model. We formulated the mass matrices 
of the neutral Higgs bosons and the charged Higgs mass in a manner independent 
of the phase convention in CP-violating case. By use of the effective potential 
including one-loop corrections of the third generation of quarks and squarks, we 
obtained constraints on the parameters of the model. Then we required that the 
neutral Higgs boson whose coupling to the Z boson is not so small must be heavier 
than 114GeV, and that the prescribed vacuum must be the absolute minimum of 
the effective potential. The latter is nontrivial in the NMSSM, in contrast to the 
MSSM, in which the electroweak vacuum is the global minimum of the potential as 
long as the symmetry breaking conditions are satisfied. We found that the vacuum 
condition leads to an upper bound on the charged Higgs boson, which is irrelevant 
in the MSSM limit but effective in our case of (N) = O(100)GeV. The allowed 
parameter regions are classified into two distinct sets: one admits a light scalar 
with very small gauge coupling, while the other contains the scalars heavier than 
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the bound. The former is realized only in the case of weak scale (N) for k ~ 0. 
In this light Higgs scenario, the light Higgs bosons are not excluded by currently 
published experimental data: the correction factors to the Yukawa coupling of the 
light bosons are not so large and those to the Z/i/i-coupling are so small that the 
light bosons have not been observed yet by experiments. Therefore the lightest 
scalar among those with large gauge coupling, which is the candidate produced in 
future collider experiments, behaves just like the Higgs boson in the standard model. 
The SM-like Higgs boson can be as heavy as 150GeV, so that the model will take 
place of the MSSM if the first observed Higgs boson is heavier than 135GeV. 

Another feature of the Higgs sector in the NMSSM is a possible CP violation at 
the tree level. We studied an explicit CP violation in the Higgs sector, which does 
not affect the neutron EDM. Such a CP violation is expected to play an important 
role in the scenario of electroweak baryogenesis. For several sets of allowed parame- 
ters in the CP-conserving case, we gradually introduced a CP phase and studied its 
effects on the masses and couplings of the Higgs bosons. As expected, the effect is 
larger for the parameter sets in the heavy Higgs region, which has larger \k\ than the 
light Higgs scenario. We, however, could not find a case that a large CP violation 
makes the lightest Higgs boson lighter than the present bound on the SM Higgs, 
while making its gauge coupling tiny not to be produced in lepton colliders. Such a 
situation has been observed in the MSSM, in which the CP violation in the squark 
sector induces a large CP violation in the Higgs sector. If such a large CP violation 
is common to all the generation of squarks, it should be constrained by the neutron 
EDM experiments. We have found that such a CP violation weakens the first-order 
EWPT in the MSSM with a light stop [13]. Thus the available parameter region 
for electroweak baryogenesis is very limited in the MSSM. As naively expected, if 
the EWPT in the NMSSM is strongly first order because of the trilinear terms in 
the Higgs potential, it will open new possibility of the baryogenesis [12]. Since the 
strong phase transition will not be caused by a light stop, it will persist in the 
presence of the CP violation in the Higgs sector. Further, the CP violation is not 
so strongly constrained by the EDM experiments. A study of the EWPT in the 
NMSSM is now in progress. 
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A The field dependent masses 

In this appendix, we summarize the field-dependent masses of the quarks and 
squarks of the third generation and those of the gauge bosons. We retain only 
the neutral components of the Higgs fields, which appear in the definition of the 
effective potential (40) and are necessary for the neutral Higgs mass-squared matrix. 
For the charged Higgs mass, we need the expressions including the charged Higgs 
fields (f)^ and but they are almost the same as those in the MSSM [13] except 
for replacement of \i with Xv n e l(f /y2. 

The quark masses are expressed in terms of the Higgs fields and vacuum expec- 
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tation values, 



ml 



2 1 j 1 2 



1 



\y b \ 2 (v 2 d + 2v d h d + h z d + ai), 



1 



m\ = \vM\ 2 = ^\yt\ 2 (vl + 2v u h u + hi + a 2 u ) 



(59) 
(60) 



where <f> d = (v d + h d + ia d )/\/2 and 0° = e ld {v u + h u + ia u )/y2. If we take the 



Higgs fields as h d = h u = h n = 0, the masses at the vacuum are obtained, 



1 



mi 



ml 



\y b \ 2 v 2 d 



m: 



m. 



\\yt\ 2 vl, 



(61) 



where the angle brackets denote vacuum expectation value. The field dependent 
masses of the gauge bosons are 



m z 



2 2 
2 _ 92 + 9l ,\ ,0|2 



#2 /| .0|2 



i.0 |2> 



(62) 



Then the gauge bosons masses are 



(m|) = m| 



52+5! 



Similarly, field dependent top- and bottom-squark masses are 



m 



'1.2 



m 



In. 



where 



x t = 7 52 - o^i 



7 52 - 7^2 



The masses of the squarks at the vacuum are 



11,2 



2, 2 ,52+51/2 2 \ , i |22 
m q + m t fl + o ( U d ~ V «) + Iftl V « 



±jM? + 2\y t \Z(P t v* + Q t vZ) 



2 i 2 52 + 5l / 2 2\ , i |2 2 
m ? + m 6fl 8 ( d ~ U > '™ d 



R q = -±=R e (\A q e^), 



1 



(63) 



m| + ml + ^±^(|^| 2 - |0°P) + 2\y t \ W 

±^(m? - m\ R + x t (\^ - |0O |2))2 + 4 | yt |2|An^ - ^^| 2 ] , (64) 
m 2 + - 9 -^(\4>° d \ 2 ~ |</>°J 2 ) + 2\y b \ 2 \$\ 2 

± " m L + - l^l 2 )) 2 + 4|5 fe | 2 |An^o _ ^0O*| 2 ] , (65) 



(66) 



±^M2 + 2|y 6 |2(n^ + Q^ ) 
where we define the following combinations of the parameters: 

i} ), (* = *,&), 



(67) 



p t = ^ I A l 2 ^ ~ ^t^n tan/?, 

P b = T^l^n - RbVn COt (3, 



I g = -^Im(\A q e^>), 
V2 

Q t = \A t \ 2 - R t v n cotp, 



Q b = \ A b \ 2 - R b v n t&nf3, 



AC 



2 2 2,- E </2 2\ 

- rriq — m~ R + —{v d — v u ), 



M 2 = mf 



ml R + ^(v 2 -v 2 u ). (68) 
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Although I q doesn't appear above equations, we defined this for later convenience 
(We use I q in Appendix B). 



B Derivatives of the squark masses 

The corrections from the squark loops to the neutral-mass matrix contain first 
derivatives of the field-dependent-squark masses and the matrices T s , T p , T SP , 
B s ', B p and B SP (50)-(56). The First derivatives of the squark masses are as 
follows, 



' dm 2 



dh 
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al+a 2 



m\ 



L v d 



2 Am 2 ' 



dm- 

*1,2 

da 



\yt\ 2 

Am 2 



IfVnP, 



' dm 2 
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2/6 1 
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2 Am?' 
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hv n p, 



(70) 



(x t M 2 + 2\y t \ 2 P t )v d 
(-x t M 2 + 2\y t \ 2 Q t )v u 




(x b M 2 + 2\y b \ 2 Q b )v d 
{-x b M 2 + 2\y b \ 2 P b )v u 
2|?/6| 2 n# 



m 



'/i 



m 



92 ' 



(q = t,b). (71) 



The detailed form of the matrices T s and B s are, 
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(72) 



bb T 



2(Am?) 
(73) 



2)2- 



These are obtained from the second derivatives of the AqV with respect to h d , h u and 
h n , and applying the tadpole conditions. The corrections to the pseudoscalar and 
scalar-pseudoscalar-mixing component of the neutral-mass matrix include following 
matrices: 



I |2 D Vn 

\y t \ Rt 

VdVu 
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and 



\yt\ 2 hvn 
(Am?) 2 



tp 




bp T . 



(75) 



References 

[1] ALEPH Collaboration, R. Barate et al, Phys. Lett. B565 (2003) 61, 
arXiv : hep-ex/0306033. 

[2] Particle Data Group Collaboration, S. Eidelman et al., Phys. Lett. B592 
(2004) 1. 

[3] Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 85 (1991) 1. 
J. R. Ellis, G. Ridolfi and F. Zwirner, Phys. Lett. B257 (1991) 83. J. R. Es- 
pinosa and R.-J. Zhang, Nucl. Phys. B586 (2000) 3, arXiv:hep-ph/0003246. 

[4] P. Fayet, Nucl. Phys. B90 (1975) 104. 

[5] J. F. Gunion and H. E. Haber, Nucl. Phys. B272 (1986) 1. 

[6] J. R. Ellis, J. F. Gunion, H. E. Haber, L. Roszkowski and F. Zwirner, Phys. 
Rev. D39 (1989) 844. 

[7] D. J. Miller, R. Nevzorov and P. M. Zerwas, Nucl. Phys. B681 (2004) 3, 
arXiv : hep-ph/0304049. 

[8] A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B553 (1999) 3, 
arXiv:hep-ph/9902371. M. Carena, J. R. Ellis, A. Pilaftsis and C. E. M. 
Wagner, Nucl. Phys. B586 (2000) 92, arXiv:hep-ph/0003180. 

[9] The OPAL Collaboration, Physics Note PN524 (2003). 

[10] U. Ellwanger, M. Rausch de Traubenberg and C. A. Savoy, Nucl. Phys. B492 
(1997) 21, arXiv:hep-ph/96 11251. G. C. Branco, F. Kruger, J. C. Romao and 
A. M. Teixeira, JEEP 07 (2001) 027, arXiv:hep-ph/0012318 A. T. Davies, 

C. D. Froggatt, and A. Usai, Phys. Lett. B517 (2001) 375. 

[11] M. Matsuda and M. Tanimoto, Phys. Rev. D52 (1995) 3100, 
arXiv:hep-ph/9504260. N. Haba, Prog. Theor. Phys. 97 (1997) 301, 
arXiv:hep-ph/9608357. S. W. Ham, J. Kim, S. K. Oh and D. Son, Phys. 
Rev. D64 (2001) 035007, arXiv:hep-ph/0104144. S. W. Ham, S. K. Oh and 

D. Son, Phys. Rev. D65 (2002) 075004. 

[12] M. Dine, P. Huet and J. R. Singleton, Nucl. Phys. B375 (1992) 625. 
M. Pietroni, Nucl. Phys. B402 (1993) 27, arXiv :hep-ph/9207227. A. T. 
Davies, C. D. Froggatt and R. G. Moorhouse, Phys. Lett. B372 (1996) 88, 
arXiv:hep-ph/9603388. S. J. Huber and M. G. Schmidt, Eur. Phys. J. CIO 
(1999) 473, arXiv:hep-ph/9809506. S. J. Huber and M. G. Schmidt, Nucl. 
Phys. B606 (2001) 183, arXiv:hep-ph/0003122. M. G. Schmidt, Nucl. Phys. 
Proc. Suppl. 101 (2001) 39, arXiv :hep-ph/0102043. 

[13] K. Funakubo, S. Tao and F. Toyoda, Prog. Theor. Phys. 109 (2003) 415, 
arXiv : hep-ph/021 1238. 

[14] B. A. Dobrescu and K. T. Matchev, JEEP 09 (2000) 031, 
arXiv : hep-ph/0008192. 

[15] S. A. Abel, S. Sarkar and P. L. White, Nucl. Phys. B454 (1995) 663, 
arXiv : hep-ph/9506359. 



20 



[16] A. Menon, D. E. Morrissey and C. E. M. Wagner, arXiv:hep-ph/0404184. 

[17] CDF and DO Working Group Members Collaboration, L. Babukha- 
dia et al, FERMILAB-PUB-03-320-E, 2003. D. Denegri et al, 
arXiv:hep-ph/01 12045. 

[18] S. Gentile, ATL-PHYS-2004-009, 2004. S. Asai et al, Eur. Phys. J. C32 Suppl. 
2 (2003) sl9, arXiv:hep-ph/0402254. K. A. Assamagan and Y. Coadou, Acta 
Phys. Polon. B33 (2002) 1347. 

[19] G. Abbiendi et al. OPAL Collaboration, Eur. Phys. J. C18, (2001) 425 
arXiv : hep-ex/0007040. 

[20] G. Abbiendi et al. OPAL Collaboration, Eur. Phys. J. C23 (2002) 397 
arXiv:hep-ex/01 11010. 

[21] V. M. Abazov et al. DO Collaboration, arXiv:hep-ex/0410062. 



21 



